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Abstract 

Background: It has been hypothesised that light skin pigmentation has arisen to ensure adequate levels of vitamin 
D as human populations moved out of Africa and into higher latitudes. Vitamin D, which is primarily obtained 
through exposure to sunlight (specifically ultraviolet radiation B (UVR-B)), has been inversely associated with several 
complex diseases. Greater sun exposure, on the other hand, is a well-known cause of skin cancer. The potential of 
UVR to be beneficial for some health outcomes but detrimental for others has prompted a public health debate on 
how to balance the positive and negative consequences of sun exposure. In this study we aimed to determine the 
validity of the evolutionary hypothesis linking lighter skin with higher vitamin D concentrations in a European 
population. Additionally, we aimed to examine the influence of pigmentation on personal behaviour towards 
sunlight exposure and the effects of this behaviour on vitamin D. 

Methods: We combined genetic variants strongly associated with skin colour, tanning or freckling to create genetic 
scores for each of these phenotypes. We examined the association of the scores with pigmentary traits, sun 
exposure and serum 25-hydroxyvitamin D (25(OH)D) levels among children of the Avon Longitudinal Study of 
Parents and Children (ALSPAC, N = 661 to 5649). 

Results: We found that fairer-skinned children, i.e. those with higher pigmentation score values, had higher levels 
of 25(OH)D (0.6 nmol/l; 95% CI 0.2, 1.0; per unit increase in skin colour score; N = 5649). These children also used 
more protection against the damaging effects of UVR. 

Conclusions: In this population taking protective measures against sunburn and skin cancer does not seem to 
remove the positive effect that having a less pigmented skin has on vitamin D production. Our findings require 
further replication as skin pigmentation showed only a small effect on circulating 25(OH)D. 
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Background 

Skin pigmentation is a complex phenotype that shows 
extensive variation between human populations. Evidence 
suggests that such variation was shaped by the action of 
natural selection at different latitudes, to prevent DNA 
damage by ultraviolet radiation (UVR) to the skin and to 
guarantee sufficient synthesis of vitamin D [1], although 
other hypotheses that explain its worldwide distribution 
have been proposed [2]. Among these are protection 
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against light-induced folate deficiency, protection from 
cold injury, prevention of heat load, camouflage and 
resistance against infectious disease [2]. The vitamin D 
hypothesis posits that as human ancestors migrated 
out of the tropics and into areas with reduced sunlight 
there was a progressive depigmentation of the skin to 
allow for UVR-B (280-320 nm) induced vitamin D 
production [1]. 

Vitamin D synthesis is highly dependent on the con- 
centration of melanin in the skin as melanin absorbs and 
scatters UVR-B, resulting in a less efficient conversion of 
7-dehydrocholesterol to previtamin D 3 [3]. Therefore, 
dark-skinned individuals will experience slower vitamin D 
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synthesis than light-skinned ones. This is especially im- 
portant at higher latitudes where the incidence and 
duration of sunlight is reduced. However, although in- 
dividuals with light skin pigmentation are more effi- 
cient at producing vitamin D, they will be affected by 
an elevated sunburn response, a lower tanning ability 
[4], and a greater susceptibility to skin cancers [5]. 

A number of studies have shown an association between 
skin lightness and 25-hydroxyvitamin D (25(OH)D) levels 
[6-10], most of them comparing individuals from different 
ethnicities. 

Other factors that reduce the cutaneous photosyn- 
thesis of vitamin D are clothing and increasing age [3]. 
Sunscreens have been found to have a similar effect 
under strictly controlled conditions but recent evidence 
suggests that their normal use may not impair vitamin D 
production [11]. 

Because solar radiation constitutes the main source of 
vitamin D, but is also a risk factor for sunburn and skin 
carcinogenesis [12], achieving a balance between gener- 
ating enough vitamin D while limiting the damage to the 
skin has become an important area of debate in public 
health (see for instance, [13]). 

Several genes underlying human skin pigmentation 
and sun sensitivity have been identified through candidate 
and genome-wide association studies (GWAS) [14-18]. As 
a result, it is possible to combine polymorphisms from 
these genes to create genetic scores that are strongly as- 
sociated with pigmentation traits, such as skin colour, 
tanning potential and freckling [19]. A score represents 
the additive effect of multiple genetic variants and is 
likely to account for a larger proportion of variance in 
the trait of interest than single genetic variants. The use 
of genetic scores instead of self- or interviewer-assessed 
pigmentation helps to avoid problems of measurement 
error and possible confounding.Using pigmentation gen- 
etic scores described earlier [19], we have examined asso- 
ciations of these scores with pigmentation phenotypes and 
sun exposure variables, and assessed their influence on 
serum 25(OH)D levels among UK children of the Avon 
Longitudinal Study of Parents and Children (ALSPAC). 
Our aim was to establish whether in this European 
population children with a lighter skin phenotype had 
higher circulating 25(OH)D as predicted by the vitamin D 
hypothesis for the evolution of skin pigmentation. Add- 
itionally, we aimed to examine the influence of pigmen- 
tation on personal behaviour towards sunlight exposure 
and whether this affects 25(OH)D levels. 

Methods 

Study population 

ALSPAC is a population-based prospective cohort study 
investigating factors that affect the health and develop- 
ment of children and their parents. The study methods 



are described in detail elsewhere [20,21] (http://www. 
bristol.ac.uk/alspac/). In brief, pregnant women living 
around Bristol, South West England, who had an expected 
date of delivery between April 1991 and December 1992 
were eligible, and of these 14,541 were enrolled in the 
study. Extensive data have been collected from the 
mothers and their offspring from pregnancy onwards by 
questionnaire, abstraction from medical notes, record link- 
age and by attendance at research clinics. Ethical approval 
for the study was obtained from the ALSPAC Law and 
Ethics Committee (IRB# 00003312) and the Local Research 
Ethics Committees (Bristol and Weston, Southmead, and 
Frenchay Health Authorities). Written informed consent 
was obtained from all participants in the study. 

Skin characteristics and sun exposure 

A 10% randomly selected sample of parents whose babies 
were born within the last 6 months of the survey were in- 
vited to bring their children to a research clinic (Children 
in Focus, CIF) at 4, 8 and 12 months old and 6-monthly 
intervals thereafter up to 61 months of age, where a num- 
ber of clinical, physiological and psychological assessments 
were completed. Skin colour was measured at 49 months 
of age with an EEL (Evans Electroselenium Ltd) spectro- 
photometer. Observers measured the skin reflectance 
on the inner surface of the arm, midway between scap- 
ula and elbow. The spectrophotometer was calibrated 
daily against a matt white ceramic tile and the interior 
of a black box to give 100% and 0% reflectances. The re- 
flectance of a yellow melamine surface was then measured 
and recorded each day. 

Detailed observations of freckles and moles were car- 
ried out at 49 months and repeated at 61 months. Be- 
fore coming to the clinic at 49 months parents were 
given an outline drawing of a child of the same age and 
sex as their own, and asked to mark on it any freckles, 
moles or other marks on the skin. At the clinic the ob- 
servers used this drawing as the focus when examining 
the child's skin. They observed the child in good light, 
naked if possible. The number and location of all freckles 
and moles visible to the naked eye was noted on body 
maps on which face/scalp, neck, trunk and pelvis, and the 
upper, mid, and lower regions of each limb were shown, 
in both front and rear view (i.e. 30 areas of the body). 
Information on having any freckles at either 49 months, 
61 months or both was combined in one exposure vari- 
able coded as 0 (no freckles at all) or 1 (at least one freckle 
at any age). Children with missing information on freckles 
at 49 months but with a valid response at 61 months were 
coded based on the latter. Children with no freckles at 
49 months and missing information at 61 months were 
coded as missing. The number of moles (acquired melano- 
cytic naevi) grouped in quartiles was used in the analysis. 
Mole count was additionally assessed during the clinic 
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visit at age 15 years. A questionnaire was sent in advance 
asking for the presence of large and small moles, as mea- 
sured by a plastic guide, to be indicated on arms and legs. 
Moles on one arm were validated on a sub-set of ~6% of 
children that attended the clinic. Total mole count was de- 
rived as the sum of the number of moles present on arms 
and legs. 

For the whole of the ALSPAC cohort questions about 
the child's previous exposure to the sun, the skin's reac- 
tion, and protective measures taken, were included in 
the questionnaires sent to parents when the child was 
around 54, 65, 69, 77, 103 and 140 months old. A com- 
bined "sunburnt" variable was created using the informa- 
tion provided as a response to the question of whether 
the child had been badly sunburnt with blisters or pain 
for at least 2 days since the previous questionnaire in- 
stance. The variable was coded 0 (not badly sunburnt from 
birth until 140 months old) or 1 (at least one occurrence 
of bad sunburn up to 140 months of age). Children who 
did not report any instance of sunburn but who were 
missing information for one or more of the questionnaires 
were coded as missing. 

Time spent outside in the sun annually either locally 
(more than 4 hours near the sea) or abroad at 5 and 
9 years old was derived from questionnaires completed 
when the child was 69 and 103 months of age. 

Additional file 1: Table SI shows measured variables, 
time points at which they were collected and numbers 
of individuals available for each time point. 

Vitamin D levels 

Serum total 25(OH)D (nmol/1) was assessed from samples 
obtained at mean age 9.9 years. If no samples were available 
from the assessment at age 9, which contained the largest 
number of children (N = 4364), samples from the 11 year 
assessment (N = 995) or the 7 year assessment (N = 1074) 
were used [22]. The age range of the children with a 25 
(OH)D measure was 85 to 163 months. 

Following collection, samples were immediately spun, 
frozen and stored at -80°C, with no further freeze-thaw 
cycles to the time of 25(OH)D analysis. 25(OH)D 3 and 
25(OH)D 2 were measured with high performance liquid 
chromatography tandem mass spectrometry using an in- 
ternal standard in a laboratory meeting the performance tar- 
get set by the Vitamin D External Quality Assessment 
Scheme (DEQAS) Advisory Panel for 25(OH)D assays. 
Inter- assay coefficients of variation for the assay were <10% 
across a working range of 2.5-624 nmol/1 for both 25(OH)D 3 
and 25(OH)D 2 . Total 25(OH)D levels were calculated as 
the sum of 25(OH)D 3 and 25(OH)D 2 [22]. 

Genotyping 

ALSPAC children were genotyped using the Illumina 
HumanHap550 quad chip (Illumina, Inc., San Diego, 



CA) by 23andMe subcontracting from the Wellcome 
Trust Sanger Institute, Cambridge, UK and the Laboratory 
Corporation of America, Burlington, NC, US. The result- 
ing raw genome-wide data were subjected to standard 
quality control methods. Individuals were excluded on the 
basis of gender mismatches; minimal (<0.325) or excessive 
heterozygosity (>0.345); disproportionate levels of indi- 
vidual missingness (>3%), cryptic relatedness measured 
as proportion of identity by descent (IBD > 0.1) and in- 
sufficient sample replication (IBD < 0.8). The remaining 
individuals were assessed for evidence of population 
stratification by multidimensional scaling analysis and 
compared with Hapmap II (release 22) European descent 
(CEU), Han Chinese (CHB), Japanese (JPT) and Yoruba 
(YRI) reference populations; all individuals with non- 
European ancestry were removed. Single nucleotide poly- 
morphisms (SNPs) with a minor allele frequency of < 1%, 
a call rate of < 95% or evidence for violations of Hardy- 
Weinberg equilibrium (p<5xl0~ 7 ) were removed. Geno- 
typic data were subsequently imputed using Markov 
Chain Haplotyping software (MACH v.1.0.16) and 
phased haplotype data from CEU individuals (HapMap 
release 22, Phase II NCBI B36, dbSNP 126) based on a 
cleaned dataset of 8,365 individuals and 464,311 auto- 
somal SNPs. 

Genotypic dosages for all pigmentation-related SNPs 
of interest, which represent the expected number of one 
of the alleles and range from 0 to 2, were used in the as- 
sociation analysis. Genotypes were checked for deviation 
from Hardy- Weinberg equilibrium using the hwsnp func- 
tion implemented in Stata (StataCorp LP, 2012, College 
Station, TX). 



Pigmentation genetic scores 

In order to minimize bias due to residual confounding 
and measurement error of self- or interviewer-reported 
pigmentation and sun exposure variables, we computed 
three genetic scores from SNPs that have previously 
been associated with the pigmentation-related pheno- 
types examined in this study: skin colour (skin colour 
score = SCS), tanning ability (tanning score = TS) and 
freckling (freckling score = FS) [14]. It should be noted 
that ALSPAC samples were not part of any of the 
GWAS that discovered the pigmentation-associated 
SNPs and therefore the scores are unlikely to be over- 
fitted in this independent sample. SNPs chosen had a 
minor allele frequency > 0.05 in the CEU population, 
and the "risk" allele had been clearly established. Scores 
were calculated by summing the dosages of all appropri- 
ate SNPs in each individual, with no missing SNP data. 
"Risk" alleles were those associated with lighter skin 
color, burning rather than tanning, and a greater likeli- 
hood of having freckles. Dosages for each SNP range 
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from 0 (when no risk alleles are present) to 2 (when the 
child carries two risk alleles). 

Polymorphisms included in each score are shown in 
Table 1. These scores have been used before in a study 
of pigmentation characteristics, vitamin D and prostate 
cancer in British men [19]. 

Population stratification 

The top 10 principal components (PCs) that reflect the 
population's genetic structure were estimated according 
to Price et al. [23] from genome-wide SNPs genotyped, 



imputed and cleaned in ALSPAC children as described 
above. All 10 PCs were included as covariates in the 
regression models to account for confounding by popula- 
tion stratification. 



Potential confounders 

We examined whether the genetic scores were associated 
with confounders. Children's date of birth and sex were 
obtained from birth records. Information on maternal 
education was collected using questionnaires completed 



Table 1 Skin colour, tanning and freckling genetic scores 


Skin colour score 












Gene 


SNP 


Chromosome 


Chromosomal location 


Major/minor alleles 


Risk allele 3 


IRF4 


rs 12203592 


6p25.3 


341321 


C/T 


T 


TYR 


rs 1042602 


11q14.3 


88551344 


C/A 


A 


HERC2 


rsl 291 3832 


1 5q 1 3.1 


26039213 


G/A 


G 


MC1R 


rsl 805007 


16q24.3 


88513618 


C/T 


T 


ASIP 


rs491 1414 


20q 11.22 


32193105 


G/r 


T 


ASIP 


rs1015362 


20q 11.22 


32202273 


G/A 


G 


Tanning score 












Gene 


SNP 


Chromosome 


Chromosomal location 


Major/minor alleles 


risk allele 13 


PPARGOB 


rs32579 


5q32 


149191041 


G/A 


G 


IRF4 


rs 12203592 


6p25.3 


341321 


CT 


T 


IRF4/EXOC2 


rs 122 10050 


6p25.3 


420489 


C/T 


T 


TYR 


rs 1393350 


11q14.3 


88650694 


G/A 


A 


CALCOCO l/HOXC 1 3 


rs7969151 


12q13.13 


52445544 


G/A 


A 


PAPOLA/VRK1 


rs 17094273 


14q32.2 


96173560 


G/A 


A 


HERC2 


rsl 291 3832 


1 5q 1 3.1 


26039213 


G/A 


G 


CPNE7/DPEPI 


rs 154659 


16q24.3 


88194838 


T/C 


C 


MC1R 


rsl 805007 


16q24.3 


88513618 


CT 


T 


DBNDD1 


rs 11648785 


16q24.3 


88612062 


C/T 


C 


ASIP 


rs491 1414 


20q 11.22 


32193105 


G/r 


T 


ASIP 


rsl 01 5362 


20q 11.22 


32202273 


G/A 


G 


PRDM15 


rs7279297 


21q22.3 


42100984 


A/G 


A 


Freckling score 












Gene 


SNP 


Chromosome 


Chromosomal location 


Major/minor alleles 


Risk allele' 


IRF4 


rs 12203592 


6p25.3 


341321 


C/T 


T 


BNC2 


rs21 53271 


9p22.2 


16854521 


A/G 


A 


TYR 


rs 1042602 


11q14.3 


88551344 


C/A 


C 


TYR 


rs 1393350 


11q14.3 


88650694 


G/A 


A 


MC1R 


rsl 805007 


16q24.3 


88513618 


C/T 


T 


ASIP 


rs491 1414 


20q 11.22 


32193105 


G/r 


T 


ASIP 


rs1015362 


20q 11.22 


32202273 


G/A 


G 


EIF6 


rs6 19865 


20q 11.22 


33331111 


G/A 


A 



a The risk allele is the allele associated with lighter skin pigmentation. 
"^The risk allele is the allele associated with higher propensity to burn. 
c The risk allele is the allele associated with a greater likelihood of having freckles. 
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by the mothers at approximately 18 and 32 weeks of 
pregnancy. 

Statistical analysis 

The main source of 25(OH)D 3 (and hence total 25(OH)D) 
is synthesis in the skin in response to UVR-B exposure, 
which means there is a strong sinusoidal pattern of 25 
(OH)D with date of blood collection. We modelled this re- 
lationship using a cosine regression model and extracted 
the residuals, using these as our main 25(OH)D measure- 
ment. These residuals remove the seasonal influence on 
25(OH)D and are a better marker of habitual 25(OH)D 
concentration for each individual, they also remove the 
noise in the 25(OH)D measure associated with the quasi- 
random sampling dates [24] . 

The association of genetic scores with pigmentation and 
sun exposure variables, and serum 25(OH)D levels, was 
assessed using one way ANOVA, t-tests, linear and logistic 
regression with adjustment for age, sex and population 
stratification. Interaction between genetic scores and sex 
on reflectance and freckling phenotypes was examined 
using a likelihood ratio test. The unique contribution of 
each genetic score to the variability of skin reflectance, 
freckling and serum 25(OH)D was calculated as the 
square of the semipartial correlation coefficients obtained 
from a model that additionally includes sex, age at clinic 
visit and population stratification principal components. 

Results 

Pigmentation genetic scores 

One SNP, rs7279297, was marginally out of Hardy- 
Weinberg equilibrium (p = 0.05). 

All scores were normally distributed. The skin colour 
score ranged from 1 to 10, the tanning score from 4 to 
19, and the freckling score from 1 to 13, a higher score 
indicating lighter skin colour, skin that burns instead of 
tanning, and a greater chance of having freckles. The 
correlation of scores after adjustment for population 
stratification was: TS vs SCS, R 2 = 0.19; TS vs FS R 2 = 
0.22; SCS vs FS R 2 = 0.01; p < 0.001 for all comparisons. 
This correlation pardy reflects the fact that some SNPs are 
present in more than one score. The correlation of scores 
with principal components representative of population 
stratification is shown in Additional file 2: Table S2. 



Potential confounders 

None of the genetic scores were associated with sex 
(Additional file 1: Table S3) or maternal education 
(Additional file 1: Table S4). The skin colour score was 
weakly associated with age of the child at completion of 
the 54 month questionnaire (mean difference per unit in- 
crease in score: 0.02 months; 95% CI 0.001, 0.05; p = 0.04), 
the other genetic scores were not associated with age. 

Sex was associated with skin reflectance, skin colour 
change after spending time in the sun, having been 
badly burnt, mole count at ~5 and 15 years of age and 
wearing protective clothing and sun block when outside 
(Additional file 1: Table S3). Some of these findings sug- 
gested that girls were slightly darker than boys (e.g. 
showing lower skin reflectance), although other results 
appeared inconsistent with this conclusion (e.g. a lower 
number of girls reported tanning easily or always tanning). 
Time spent in the sun at ~5 and ~9 years old, whether 
in Britain or overseas, did not differ between girls and 
boys (Additional file 1: Table S3). Sex was associated 
with serum 25(OH)D, with girls having on average 
lower levels than boys (mean difference: -3.41 nmol/1; 
95% CI -6.34, -0.48; p = 0.02, adjusted for age at blood 
draw). There was a similar difference between sexes after 
additional adjustment for skin reflectance and sun block 
use (mean difference: -3.53 nmol/1; 95% CI -6.52, -0.55; 
p = 0.02). 

Skin colour change after being in and out of the sun 
for a few days, experiencing a bad burn, wearing protect- 
ive clothing and sun block, and time spent outside lo- 
cally or abroad were associated with maternal education, 
with a tendency for children of better educated mothers 
to burn easily (although having a bad burn less often), 
frequently use sun protection and spend more time in 
the sun. In addition, circulating 25(OH)D was marginally 
lower for children of mothers with a high educational 
achievement (Additional file 1: Table S4). 

Pigmentation genetic scores and sun exposure variables 

Pigmentation scores were robustly associated with skin 
colour (as measured by skin reflectance, Table 2), tan- 
ning ability, skin reaction, freckling and mole count 
(Table 3), in the directions expected. Moreover, scores 
were associated with patterns of behaviour that may re- 
flect having a certain skin complexion, such as wearing 



Table 2 Pigmentation scores and skin reflectance 

Pigmentation scores Effect 3 (unadjusted) 95% CI p-value R 2 R 2 (%) F Effect 3 (adjusted b ) 95% CI p-value 

Skin colour score 0.41 0.25,0.57 8.1x10~ 7 0.035 3.46 24.8 0.41 0.25,0.57 7.6x10" 7 

Tanning score 0.19 0.10,0.28 5.6x1 0" 5 0.023 2.32 16.4 0.18 0.09,0.27 9.3x1 0" 5 

Freckling score 0.25 0.13,0.36 1.9x10~ 5 0.026 2.62 18.6 0.22 0.11,0.34 1.0x10~ 4 



a Mean difference in skin reflectance per unit increase in pigmentation score. 
b Adjusted for sex, age at clinic visit and principal components. 
N = 693. 
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Table 3 Distribution of reported or observed pigmentation characteristics and 25(OH)D levels by pigmentation genetic 
scores 



Mean ± SD of pigmentation score 
Skin colour score Tanning score 



Freckling score 



Mean ± SD 
25(OH)D (nmol/l) 



Skin colour change 

Always burns never tans 
Burns easily rarely tans 
Doesn't change 
Tans easily rarely bums 
Always tans never burns 
p-value 

Child badly burnt 

No 
Yes 

p-value 
Any freckles 

No 
Yes 

p-value 
Mole count 
49 months 

Q1 (0-6) 
Q2 (7-9) 
Q3 (10-14) 
Q4 (15-57) 
p-value 
61 months 
Q1 (0-10) 
Q2 (11-16) 
Q3 (17-25) 
Q4 (26-86) 
p-value 
186 months 
Q1 (0-2) 
Q2 (3-9) 
Q3 (10-22) 
Q4 (23-153) 
p-value 



6.29 ± 1 .25 

5.69 ±1.27 
5.41 ±1.18 
4.89 ±1.17 
4.57 ±1.11 
<0.0001 

5.00 ±1.21 
5.40 ± 1 .29 
1.5x10~ 13 

4.70 ±1.14 
5.45 ± 1 .23 
2.0x1 0" 1 6 



4.84 ±1.17 
4.94 ±1.18 
5.03 ± 1 .25 
5.37 ±1.22 
0.0002 

4.83 ± 1 .22 

4.84 ± 1 .24 
5.05 ±1.1 7 
5.30 ±1.25 
0.001 

5.05 ± 1 .22 
4.92 ± 1 .20 
5.00 ±1.21 
5.21 ±1.18 
0.0001 



1 2.1 6 ± 1.97 
11. 76 ±2.1 9 
11.1 9 ±2.02 
10.35 ±2.00 
9.92 ± 1 .95 
<0.0001 

10.54 ±2.08 
11. 27 ±2.1 7 
2.2x1 0" 15 

1 0.06 ± 1 .96 
11.41 ±2.18 
2.1x10" 17 



1 0.27 ± 1.90 
10.30 ±2.04 
10.75 ±2.22 
11. 24 ±2.27 

I. 89x10~ 5 

10.07 ±2.06 
10.43 ±2.11 
10.68 ±2.1 6 

II. 21 ±2.10 
3.50x10"' 

10.56 ±2.1 5 
10.50 ± 2.14 
10.60 ±2.08 
10.83 ±2.04 
0.02 



7.1 6 ± 1.54 
6.82 ± 1.81 

6.41 ±1.66 
5.85 ±1.61 
5.56 ± 1.54 
<0.0001 

5.96 ± 1 .68 

6.42 ± 1 .70 
6.4x10" 10 

5.52 ± 1.57 
6.79 ± 1 .80 
2.3x1 0" 22 



5.80 ± 1 .60 
5.66 ± 1 .64 
6.08 ± 1 .95 
6.49 ± 1 .80 

3.1 6x1 0" 5 

5.55 ±1.75 
5.90 ± 1 .67 
6.23 ± 1 .82 
6.43 ± 1 .73 
4.26x10" 1 

6.05 ± 1 .76 
6.03 ± 1 .72 
6.01 ±1.67 
6.03 ± 1 .65 
0.99 



■-1/ 
5--1/ 

235 
2209 
924 
3962 

3501 
613 
4114 

432 
297 
729 



227 
144 



176 
733 



191 
1 74 



166 
711 



691 
710 
665 
2754 



62.8 ±14.6 
63.4 ± 1 6.7 

64.9 ± 1 8.4 
63.2 ± 19.9 

63.6 ±20.1 
0.81/0.53 c 

62.4 ± 1 8.7 

63.7 ± 19.3 
0.1 3/0.1 T 

60.8 ±19.1 
63.7 ±19.2 
0.05/0.04 c 



61.3 ± 16.6 
60.2 ± 1 5.4 
61.5 ±25.2 
62.8 ± 1 6.7 
0.71/0.74 c 

60.4 ± 1 7.6 

60.5 ±21.7 
61.4 ±15.1 

65.2 ±21.5 
0.08/0.09 c 

60.6 ± 1 7.0 

61.4 ±19.1 

62.3 ±18.5 

63.5 ±18.5 
0.02/0.01 c 



45 

518 

225 

2098 

875 

3761 

3465 

594 

4059 

391 
270 
661 



193 
123 
166 
169 
651 

164 
162 
163 
160 
649 

704 
717 
729 
656 
2806 



a N = number of individuals included in the analysis of pigmentation scores. 
b N = number of individuals included in the analysis of 25(OH)D levels. 
c Adjusted for sex and age at blood draw. 

protective clothing or applying sun block (Additional 
file 1: Table S5). The effect of the freckling score was 
weaker with respect to the behavioural variables com- 
pared to the effects of the skin colour and tanning 
scores. Children with higher scores (i.e. lighter skin 
colour, who tend to burn and have freckles) were more 
likely to always be covered or use sun block when out in 



the sun than children with lower scores. Adjustment for 
population stratification did not alter the results appre- 
ciably (data available on request). 

Except for the tanning score and spending time in the 
sun when abroad at 9 years of age (mean difference in tan- 
ning score for children who spent time in the sun -0.35; 
95% CI -0.61, -0.08; p = 0.01), there were no associations 
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between the pigmentation scores and whether children 
spent time in the sun (data available on request). 

There was no clear evidence of an interaction between 
the skin colour score and sex on skin reflectance (p for 
interaction = 0.12), but we found some evidence of an 
interaction between the freckling score and sex on having 
any freckles (p for interaction = 0.02). The association of 
the genetic score with the pigmentation phenotype was 
stronger in girls for all scores and traits, with the excep- 
tion of the effect of the freckling score on reflectance 
which was considerably larger in boys (Additional file 1: 
Tables S6 and S7). The unique contribution of each gen- 
etic score to total variance in skin reflectance and having 
any freckles was between 1.5 and 4 times greater in girls 
than in boys (Additional file 1: Tables S6 and S7). 

Sun exposure variables and vitamin D levels 

Children who were reported to never use sun block 
when out in the sun exhibited lower 25(OH)D levels 
compared to children who always did so (Additional file 1: 
Table S5). Children with freckles, and with a greater num- 
ber of moles at ages 5 and 15 (Table 3) had higher serum 
25(OH)D as did those who spent time in the sun abroad 
(Additional file 1: Table S8). There was no association 
of skin colour change after being in the sun, having ex- 
perienced a bad sunburn, mole count at age 4, wearing 
clothing when outside, and sun block factor used, with 
25(OH)D concentration (Table 3 and Additional file 1: 
Table S5). Likewise, skin reflectance was not observa- 
tionally associated with 25(OH)D, (mean difference per 
unit increase in reflectance: 0.12 nmol/1; 95% CI -0.45, 
0.68; p = 0.68; adjusted for sex and age at blood draw). 

Pigmentation genetic scores and vitamin D levels 

The pigmentation genetic scores were associated with 
serum 25(OH)D, i.e. the lighter the skin (as indicated by 
the genetic score) the higher the circulating 25(OH)D 
levels (Table 4). Additional adjustment of each score for 
the other two showed an association of the skin colour 
score only (mean difference per unit increase in skin 
colour score: 0.5 nmol/1; 95% CI 0.1, 1.0; p = 0.03). Scores 
explained a small proportion of the variability in 25(OH)D 
levels (R 2 < 0.002). 



Table 4 Pigmentation scores and 25(OH)D levels 



Pigmentation scores 


25(OH)D (nmol/l) a 


95% CI 


p-value 


R 2b 


Skin colour score 


0.6 


0.2, 1 .0 


0.003 


0.0015 


Tanning score 


0.3 


0.1, 0.6 


0.003 


0.0014 


Freckling score 


0.4 


0.1, 0.7 


0.01 


0.0012 



a Mean difference in 25(OH)D levels per unit increase in pigmentation score, 
adjusted for sex, age at blood draw and principal components. 
b Semipartial correlation coefficient squared, indicates the unique contribution 
of the genetic score to trait variability. 
N = 5649. 



Discussion 

Our results show that in this European population fairer 
skin was linked to higher concentrations of 25(OH)D, in 
agreement with the vitamin D hypothesis for the evolu- 
tion of skin pigmentation. We used genetic scores that 
are robustly associated with pigmentation traits and sun 
exposure measurements in order to reduce bias due to 
self-report or clinical assessment. Whilst age, sex and 
maternal education were all associated with one or more 
reported or assessed pigmentation phenotypes and also 
associated with 25(OH)D, they were mostly unassociated 
with our genetic scores, supporting the use of genetic 
scores as unbiased and unconfounded proxies for pig- 
mentary features. In cohorts with genotypic data avail- 
able but no sun exposure measurements or low quality 
sun exposure information we could use genetic scores to 
assess the relationship between sun exposure and an 
outcome of interest. 

In our study population higher scores corresponded 
to higher skin reflectance (i.e. lighter skin colour), a 
greater likelihood of having freckles, a higher mole 
count, a skin that always burns and never tans, having 
experienced at least one bad sunburn, always covering 
the skin with clothing or sun block when out in the 
sun, and using sun lotion with a high sun protection 
factor (SPF). The skin colour score, for instance, ex- 
plained ~3.5% of the variability in skin reflectance 
among ALSPAC children. 

The association of the genetic scores with sun expos- 
ure phenotypes confirms previous results obtained in 
British men who participated in the Prostate Testing for 
cancer and Treatment (ProtecT) study, where the same 
pigmentation scores were used [19]. 

The pigmentation scores, in particular the skin colour 
score, were associated with 25(OH)D levels, with higher 
scores predicting higher serum 25(OH)D, despite being 
associated with behaviours that would reduce the impact 
of UVR exposure on vitamin D synthesis. However, the 
magnitude of the effect was small. Although darker chil- 
dren (i.e. those with lower pigmentation scores) had 
lower 25(OH)D levels on average, only 0.5% of children 
were vitamin D deficient in this population (serum 25 
(OH)D < 25 nmol/1 [25]). These results do not concur 
with our findings in the ProtecT study which showed 
that British men with lighter skins -as suggested by their 
genetic scores- had lower circulating 25(OH)D [19], or 
with a study of British women in which lighter skin pig- 
mentation (assessed by nurses using the Fitzpatrick scale) 
was also associated with lower levels of 25(OH)D [26]. 
These individuals possibly used more abundant sun pro- 
tection or restricted their exposure to the sun. Children 
not wearing enough or effective sun protection and/or 
spending more time outside than adults may explain the 
reverse association. 
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There were clear differences by sex in reported sun 
exposure patterns as well as in pigmentation traits that 
suggested that girls were to some extent darker than 
boys. Yet, although apparently more protected from 
sunburn, girls reported more consistently applying sun 
block (or having sun block applied on them) when out- 
doors than boys. Time spent outside in the sun was 
similar for boys and girls. Girls have been found to have 
darker skin than boys in pre-adolescence, with boys be- 
coming darker than girls around puberty [2,27]. In our 
study there were more girls than boys with a low mole 
count at 4-5 years of age but this pattern appeared to 
have reversed by 15 years. In adults, across a range of 
countries and ethnicities, men tend to be more pigmen- 
ted than women [2,28], however, lighter-skinned men 
were identified in a study carried out in four European 
countries[18], in a sample of European Americans [29], 
and amongst populations in the Netherlands, Belgium 
and the Basque Country [20,21]. 

On the other hand, mean pigmentation genetic scores 
did not vary between sexes, suggesting that the observed 
gender differences in pigmentation phenotypes may be 
explained by environmental, hormonal or other genetic 
causes. For instance, it is well-known that cutaneous pig- 
mentation is modulated by growth factors, cytokines and 
endocrine hormones, which stimulate melanocyte prolif- 
eration anchor melanogenesis [30,31]. 

The association of the pigmentation scores with skin 
reflectance and having freckles was for the most part 
stronger in girls (with a higher reflectance and likelihood 
of freckling as the pigmentation score increases towards 
the lighter phenotype), accounting for about twice the 
phenotypic variance explained by boys' scores. A possible 
reason for this finding may be sexual selection by men for 
lighter-skinned women during human evolution [27,32]. 
In fact, a few of the SNPs included in the pigmentation 
scores are located in genes that have been targets of selec- 
tion (i.e. MC1R, HERC2, TYR, ASIP) [33-35]. However, be- 
cause we concurrently found lower reflectance values in 
girls, other factors could be influencing skin colour in the 
opposite direction in the pre-pubertal stages. A recent 
study on mate preferences among UK twins revealed a pa- 
tent sexually dimorphic preference for skin colour but, un- 
expectedly, women preferred fair-skinned men whilst men 
showed a preference for women with olive skin. Although 
a familial environmental effect was uncovered, no genetic 
effects on mate preference for skin colour were evident 
in this study [36]. These sex differences deserve to be 
examined more carefully in the future as they may have 
an impact on the onset and progression of melanoma 
[37]. Alternatively, given our small sample size and the 
lack of strong evidence for statistical heterogeneity be- 
tween girls and boys, the observed sex disparities may 
simply be due to chance. 



Differences in sun exposure-related variables were also 
detected with respect to maternal educational achieve- 
ment, such that children of better educated mothers were 
more likely to burn rather than tan, wear protective cloth- 
ing and use sun block, as well as were more often exposed 
to the sun in Britain or abroad. These children were also 
found to have slightly lower 25(OH)D levels. These find- 
ings possibly reflect reporting bias since educated mothers 
are usually more aware of the risks and benefits posed by 
the sun, and can also afford sun protection and summer 
holidays. Similar results with respect to educational level 
were reported in an observational study that assessed the 
relationship between sun exposure habits, sun protection 
behaviour and readiness to increase sun protection in an 
adult primary health care population in Sweden. Subjects 
of lower educational attainment reported being sunburnt 
less frequently, used less sun protection and revealed a 
lower propensity to increase sunscreen use [38]. 

Limitations 

Some of the pigmentation and sun exposure variables 
were measured only once in children between birth and 
age 17 (latest data release on ALSPAC children at the time 
of analysis) and therefore, it was not possible to assess 
their variation with age or to examine the contribution of 
pigmentation genes to these exposures at different ages. 
Other variables were measured at time points that were 
close together so were not very informative about age- 
related variation either. Time spent outside was only ana- 
lysed based on having holidays near the sea or abroad and 
not with basis on the children's regular outdoor experi- 
ence. Given that in the 90s, when ALSPAC children were 
growing up, there was less awareness of UVR exposure 
risks and the industry of sun protection was less developed 
-hence the low sun block SPFs available- parents' behav- 
iour towards their offspring sunlight exposure may have 
been different from that of present-day parents of young 
children [39] . Therefore, findings from this study may not 
be fully generalisable to the present day. 

Because some SNPs were associated with more than 
one pigmentation trait, they were included in more than 
one score and therefore, scores were correlated. Scores 
explained some of the variation in the non-specific traits, 
but the strongest association was seen with the specific 
phenotypes. Since we were mostly interested in assessing 
the influence of a lighter skin phenotype overall on sun 
exposure measurements and 25(OH)D levels, and not to 
unequivocally establish the effects of each specific pigmen- 
tation trait, the correlation between scores does not bias 
our conclusions. However, we acknowledge that to under- 
stand the role of specific pigmentation traits (i.e. skin 
colour, tanning and freckling) on the efficiency of vitamin 
D synthesis other approaches would be needed [40,41]. 
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We had limited statistical power to examine the as- 
sociation of some of the sun exposure variables with 
the genetic scores and with 25(OH)D levels, especially 
with respect to the interaction analyses, consequently 
the evaluation of these associations in larger studies is 
essential. 

Conclusions 

We have found that children with a lighter skin phenotype 
(light skin colour, burning, and freckles, as suggested by 
their genetic scores) had higher circulating 25(OH)D than 
their peers, in spite of the fact that they were more likely 
to wear sun protective clothing or sun block. These find- 
ings concur with the vitamin D hypothesis for the evolu- 
tion of skin pigmentation, and indicate that the benefits of 
having a less pigmented skin at a northern latitude do not 
seem to be eradicated by the use of sunscreen or clothing 
in UK children. However, we did not have enough infor- 
mation to establish whether the sun avoidance behaviour 
of ALSPAC children was indeed protective from UVR 
damage. If greater protection is deemed necessary in 
children to prevent sunburn and skin cancer, its effects 
on 25(OH)D levels may become more relevant. In con- 
trast, as reported in an earlier study, UK white men with 
lighter skin showed lower concentrations of 25(OH)D 
[19]. Our results suggest that only a small proportion of 
the variation in 25(OH)D is explained by the pigmenta- 
tion genetic scores and thus require replication in other 
studies of children before we can conclude that they are 
robust. 

Lastly, we have also shown that personal behaviour to- 
wards the sun is affected by an individual's pigmentary 
features, for instance, children of lighter complexion 
were more likely to always cover their skin with clothing, 
always wear sun block, and to apply sun block with an 
SPF 25 or higher. 
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